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Abstract: We present precise asymptotics of the odometer function for the 
internal Diffusion Limited Aggregation model. These results provide a better 
understanding of this function whose importance was demonstrated by Levine 
and Peres [3j. We derive a different proof of a time-scale result by Lawler, 
Bramson and Griffeath [2]. 
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1 Introduction 

The internal Diffusion Limited Aggregation model was first introduced by 
Diaconis and Fulton in [1] gives a protocol for building a random set recur- 
sively. At each step, the first vertex visited outside the cluster by a simple 
random walk started at the origin is added to the cluster. The resulting 
limit shape is the euclidian ball, as proved in 1992 by Lawler, Bramson and 
Griffeath in [2]. 

More recently, however, Levine and Peres [3], [4| have shown that this 
model is related to the rotor-router and divisible sandpile models. In the 
former, random walkers are replaced by eulerian walkers. In the latter, each 
vertex can hold 1 unit of mass, and the excess is divided equally among 
its neighbors when the vertex topples. Thus an initial mass at the origin 
becomes a stable shape after a suitable infinite series of topplings. 

For each of these three models, one can define an odometer function, 
which will be the total number of times a walker passes through a given 
point (counting multiple passages of the same walker) for the internal DLA 
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and rotor-router model, and the total mass emitted from a given point 
in the construction of the cluster for the divisible sandpile model. This 
function plays a comparable role in all three models and turns out to be 
instrumental in their relation. While the limiting shape of these models is 
known to be the Euclidian ball, the behavior of their odometer functions 
in the particular case where all the mass is started at the origin remains to 
be studied. 

This paper provides a closer look at the odometer function in the case of 
the internal DLA model, with an almost sure convergence of the normalized 
functions and asymptotics of this function near the origin. These results 
provide in turn a new proof of the time scale of the cluster introduced in 

m 



2 Definitions, main results 

Let (S^jgpj a sequence of independent simple random walks on Z d , and 

let us define the cluster A(n) and stopping times ((Jk)km recursively in the 
following way: 

(T = 0, 

A(0) = {0} = {SVo)}, and for all j > 0, 

Vj>0, a, = w£{t>0:Sl(t)?A(j-l)}, 

A(j) = Atf-lJU^fo)}. 

Let ||.|| denote the Euclidean norm on M d , B r the Euclidean ball of 
radius r of Z d , and Ud the volume of the unit ball of M. d . We will consider the 
cluster A{ijjdn d ), which has the same volume as B n . Lawler, Bramson and 
Griffeath proved in |2J that the normalized cluster -A( \ud,n d \ ) converges 
to the Euclidean unit ball with probability one. Lawler then improved this 
result, defining the inner and outer errors as follows: 

5j(n) = n — inf | \z\ |, 

z^A( [w d n d \ ) 

So(n) = sup \\z\ \ — n. 

zeA([w d n d \) 

to get the following bounds, with probability 1: 
lim ~~tjtt\ = °' and 

lim fofr) = 0. (2.1) 

An important part in the study of this model is played by Green's 
functions for the simple random walk, which we introduce now, as defined 
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in [5]. Let S x be a simple random walk started from a point x G Z d , and 
be the time at which S exits B n , that is : 

= mf{t : S t i B n }. 

Then we define for all y G 7L d 

Gn(x,y) = ^ 1 S x {t)=y 
t=0 

This is zero if either x or y lies outside B n . For d > 3 we define, for all 

x,y £ Z d , 

oo 

G(X,3/) = ls«(t)=y. 
i=0 

Let z G M d be a non-zero point in the open unit ball. Let us write 
z = (xi, • • • 

Let us define z n := ( [^XiJ , • • • [^dj ) • 

The odometer function m„(z) at rank n measures the number of walkers 
passing through z n in the process of building the cluster A([uj d n n \). It is 
defined as follows: 

u n {z) := Yl 1 ^{t)=z n 

i=0 t=0 

Our results are as follows: 

Theorem 2.1 For all non-zero points of the open unit ball z, when n goes 
to infinity, 



U ni Z ) II.M2 



n 2 



z\\ — 1 — 2 ln( | \z\ |) almost surely if d = 2, 



^ I|Z|12 + (rf-2)^|| rf - 2 ~ d^2 dm0St SUTdy ifd - 3 - ( 2 - 2 ) 



Remark: In [2], the authors estimate the time it takes to build a clus- 
ter of radius n, that is to say the total number of steps done by random 
walks during the construction of the cluster. The authors count these steps 
by estimating the number of steps for each given random walk. Our con- 
vergence result [27X1 allows us to take a different perspective on the problem, 
and count the total number of steps as the sum over all points in the cluster 
of the number of steps through this point. 

This leads to a different proof of the result of Lawler, Bramson and 
Griffeath on the time scaling of the cluster, which is presented in section El 
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The functions in 12.21 vanish when \\z\\ tends to 1, and tend to infinity 
when | \z\ | tends to 0. To understand this behavior around the origin, we can 
use a different scaling. Consider a sequence of non-zero points y n such that 

converges to as n tends to infinity. If we now consider the following 
sequence of values takes by the n-th odometer at point — : 

i=0 t=0 

We get the following result: 

Theorem 2.2 Let y n be a sequence of non-zero points such that ^ con- 
verges to 0. 

• If y n converges to a point a G 7h d , d > 3, when n tends to infinity, 

u n (— ) ~ u}dn d G(0, a) almost surely 
n 

• If \ \y n \\ tends to infinity, when n goes to infinity, 

almost surely, if d = 2, 

almost surely, if d> 3. 

Remark: The question of the asymptotics of the odometer function 
near the boundary remains open, and is presumably linked to the difficult 
problem of the fluctuations of the cluster around its limiting shape. 

Let us comment further on Theorem 12.11 and give a heuristic of it when 
d > 3 based on the work of Levine and Peres ([3]) (the two-dimensional 
case is similar if a little more technical). In the case of the divisible sandpile 
model, each site of Z d contains a continuous amount of mass. A site with 
an amount greater than 1 can topple, that is to say keep mass 1 and 
distribute the rest equally between its neighbors. With any sequence of 
topplings that topples each full site infinitely often, the mass approaches a 
limiting distribution, which does not depend on the sequence. 

The odometer function u for this model is defined as the total mass 
emitted from a given point. Since each neighbor y of a given point x 
divides mass equally between its 2d neighbors, the total mass received 
by x is -^Yliy^x u ^y) ■ If we define the discrete Laplacian A as Af(x) = 
55 ^y~x f(v) ~f( x ), we S et that 

Au(x) = u(x) — cr(x), 




~ 4 In 

n 2 

M v it) 2 f\\ Vn : 

n 2 d — 2 \ n 
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where a and v are the initial and final amounts of mass at x, respectively. 

In our case, the initial mass u>dn d is concentrated at the origin, and the 
final mass is 1 in each fully occupied site of the cluster. Hence: 

n — 1 if x = 0, 
Au n (x) — ^ — 1 if x is a non-zero point in the cluster, 
outside the cluster. 

To solve this equation, we first introduce a function j n that has discrete 
Laplacian a n — 1, in our case: 

ln {z) = -\\z\\ 2 - u d n d G{0, z). 

Then the odometer function u n is given by lemma 3.2 of [4] as the 
difference between the least superharmonic majorant of 7„ and 7„: 

$ n 7»u 

where s n = inf{/(x) such that / is superharmonic and / > 7„}. 

Recalling the notations z = (xi, ■ • • ,Xd) and z n = ([nxi\, ■ ■ ■ [nxd\), 
we get: 

1 2 

— 7n(*n) -> 1(Z) := ~\\Z\\ 2 - -T^WA?^ 

n z a — 2 

The function 7 is radial, and its particular shape makes it easy to determine 
its superharmonic majorant s, as one can see on Figure 1: 



j(z) if \\z\\ > 1. 



Since 7 n converges to 7, s n should converge to s. This would mean that 
u n converges to s — 7. This method is in fact made rigorous in [3], with 
the noticeable difference that in the case considered the starting mass has 
a bounded density, as opposed to our case where all of the mass is started 
at the origin. The authors then prove that the odometer function for the 
divisible sandpile model is the expected value of the odometer function 
for the internal DLA model, which gives a heuristic proof of the following 
result: 

E (Mf) ^ M 2_ 1 _ 2ln(M) if d=2, 



\ n2 / 

E ( ^f] HMI 2 + n t nd ■ — -7-^ — if d > 3. 

\ n 2 j 11 11 (d - 2)\\z\\ d - 2 d-2 

This analytic method could probably be rendered rigorous, but we prefer a 
more probabilistic approach based on the convergence of the cluster, which 
gives us an almost sure convergence result. 
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3 Proof of Theorem 2.1 



The main idea of this proof is to use the fact that the cluster is very 
similar to a disc when n is large enough to get two inequalities framing 
u n {z) between random valuations which do not depend on the shape of the 
cluster. 

Lawler's theorem proves that the following set inequality holds for n 
large enough with probability one: 

23 n _ n i/3(i nn )2 C A(\uj d n d \) C £> n+n i/3(i nn )4, 

Let us define: 

rjj(n) = n 1//3 (km) 2 
Vo{ n ) — n l ^ihin) A 

Let us define the stopping time £ J k as the time at which S- 7 leaves a ball 
As a consequence of the set inequality, we get the following inequality 

on stopping times for k large enough, for all j such that \yjd{k — \) d \ < 

j < [w^J ' w ith probability one: 
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We get the following framing of u n (z): 

n-1 L w d( fc + 1 ) d J 5fc-»j/(fe) n [u d {k+l) d \ ?fe +7)o(fc) 

Y Y Y 1 sKt)=z n < u n (z) <Y Y Y M*)-*. 

This inequality holds for n large enough with probability one because 
1 1 z n \\ — n\\z\\ is smaller than a constant. This ensures that the inequality 
(13. IB will start to hold as soon as k becomes greater than for instance. 
All the other terms of the sum are asymptotically zero. 

We will now show that the left and right bounds on u n {z), which we 
will call u n (z)~ and u n (z) + , once normalized, converge almost surely to the 
same function. In order to do this, we take a look at the following family 
of random variables: 

Xi, n ( Z ) = Y 1 SHt)=z n 
t=l 

A direct application of the Markov property shows that for i > 1, there 
exists po and p depending on k, n and z such that: 

nxi n {z) = i)=poP i -\i-p) 

Since we know that ^(xini 2 ) = 0) = P(r 2ji > where r Zn is the 
reaching time of z n , and that ^(Xk n ( z )) = Gk(0,z n ), where G n is Green's 
stopped function, a simple computation shows that: 

Po = P(6<7v„) 

n < r Zn ) 



p 



Cfc(0, z n ) 



This determines the parameter of the geometric distribution x{ n ( z ) 
follows. 

We will need the following lemma which gives an estimate of G n (0, z): 
Lemma 3.1 If z G B n , z ^ 0, we have 

2 n 1 1 

GJO, z) = — In — + O(-n-n-) + O(-) ifd = 2, and 

u>2 \\z\\ \\z\\ n 

G n (0,z) = ^-(Ikir^-^ + Odl^l 1 ^), 

where the O are uniform on the unspecified variables. 

This lemma is due to Lawler and can be found in section 1.6 of [5]. 
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We will now prove that ^m„(z) converges almost surely to the limit of 
its mean value. We first compute the mean value of u n {z)~ , then we will 
bound its variance. 

Consider the left side of (J3j) : 

Un(z)~ = E E 1 S^(t)=z n 

n-l M fe+1 ) 2 J 

Un{z)- = E **-„,(*),«(*) 

k=0 j=[uj d k 2 i 

„_1 L^2(fe+l) 2 J 

E(u n (z)~) = E G k - Vl {k)(0,Zn) 

n-l 

= ([Mk + l) 2 \-[^k 2 \)G k ^ l{k) (0,z n ) 

k=[\\z\\n\ 

(3.2) 

Let us consider separately the case d — 2. In that case, we have: 

%W1 = E (^(A: + l) 2 J-M 2 J)-fln^^ + 0(-)') 

E(u n (z)-) = E ( 4Hn ^) +°( ri5/3 ) 

^EK( 2 )-) = 1 V ^-ln-^^+OCn- 1 ^) 
rr n \ n n\\z\ \ I 

fe=L||z||nJ V 11 117 

(3.3) 

We recognize a Riemann sum, thus : 

±E(u n (z)-) -> I' 4tln(-±-)dt=\\z\\ 2 -l-21n(\\z\\) 

n J\\z\\ \\Z\\ 

When d > 3, we get instead the following estimate for E(u n (z)~) : 

n-l 

(\,,,Jk 4- 1 V* I - I 

d — 2 



E ( M fc + !) d J - Ml )tt9 77((»l 1^1 l) 2 ^-(^-^(^)) 2 ^)+o(n 1 - d ) 



fc=L||*||nJ 



8 



EM*)") = I (^^- 1 ((n||^||) 2 ^-(A;) 2 ^)')+0(n 8 / 3 - rf ) 

U k=l\\z\\n\ ^ ' 
I — I ll-ILI \ ' 



k= [I \z\ \n\ 



Z" 1 2d 

./lull ^ ~ 



EM*)~) - / ^(^||z||'-*-f)df=|MI 2 + 



n 2 v ™ w ' 7n 2| |rf-2 v 11 11 ' {d-2)\\z\\ d - 2 d-2 

The same arguments show that ^2K(u n (z) + ) has the same limit as its 
counterpart. 

We will now give a bound on the variance of these variables. Let us 
first state the following lemma: 

Lemma 3.2 The random valuations X kn ( z ) have the following variance, 
for all j, k,n G N, z G 7L d provided that \ \z\\n < k: 

Var(xi, n (z)) =(l- ^^l^ ) G k (0, z n )G k (z n , z n ) ( 2 + 1 ) -G k (0, z n f 

V G k (\J,z n ) J \ G k (z n ,z n ) - 1 J 

The proof of this lemma is a straightforward calculation which relies 
only on the knowledge of the law of x{ n ( z )- 

We will need bounds for G k (z n ,z n ). Provided that \\z\\n > k, z n is in 
the ball of radius k centered at the origin, and it has at least one neighbor 
in this same ball, which proves that G k (z n , z n ) > 1 + -^p- 

When d > 3, we can use the simple upper bound given by G k (z n , z n ) < 
G(0,0) < co. 

However, when d — 2, we will use the following bound: G k (z n ,z n ) < 
G n+k (0, 0). This bound comes from the fact that a random walk started in 
z n will exit the ball of radius k centered in before it exits the ball of radius 
n + k centered in z n . When k < n, we will even use G k (z n , z n ) < G2n(0, 0). 

We are now ready to bound the variance of u n {z)~ . 

V-p^M = ^E E Var(xU«») 

~ E E G 'fe-r, / (fe)(0, ^n)G' fc _^( fc )(z n , z n ) 2 + — 7- : — - 

n ^ i= ^j V ^fc- W (fc)K,^»;-v 

When d > 3, naming .fQ a suited constant depending only on the 
dimension d, we get: 

^ 7 k=0 
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This formula is the same as in the estimation of K(u n (z) ) and yields: 

Vai . ,Mz)_\ < g(\\z\\,d) 



n 2 J n 2 



where g is a function of ||z|| and d only. 

When d = 2, taking a suited constant K 2 gives us: 

/u^^j < _Lg^ Gfc(0 ,^)ln( 

^ ' k=0 



n) 2 , 



which yields: 



u n (z 



Var < 



Hn) 2 h(\\z\\) 



n 2 ) n 2 



where h is a function of ||;z||. 

In both cases, the sum Yl Var j is finite, which means 



M z ) \ ; D fi™^ ^u;^ ™„<™ 

n 2 

converges almost surely to lim n ^oo SH^y) ) _ The same set of arguments can 
be used to prove that Un ^ converges almost surely to Hindoo SHaMU _ 

Since lim^JMll = i im _ and u -( 2 ) < ^ < u +( z ) 

for all n large enough, we have proved our result. 



4 Proof of Theorem 2.2 



The proof of the first assertion relies only on the fact that the following 
framing holds for n large enough: 

Y Y Y 1 ^{t)=y n < Y Y 1 sKt)=yn 

fc=LlnnJ j=\w d k d \ t=1 i=\u; d (\nn) d \ 4=0 

n [uj d (k+l) d \ Cfc+, lo (fc) 

- Y Y Y 1 SKt)=y n 

fc=l+LlnnJ j=\u d k d \ t=1 

This is true because the equation (13.11) holds as soon as k is large enough, 
which happens eventually since k > Inn. 

The difference to the relevant quantity can be bounded using this in- 
equality which holds for n large enough: 

Y Y lsi ®=v™ - Y Y 1 &{t)=yn 

i=0 t=0 i=0 t=l 
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It is a consequence of the calculations in the proof of Theorem 12. II that 
the right side tends to zero almost surely once divided by n 2 . 

Just like in the proof of Theorem 12.11 we estimate the expected value 
of our lower bound: 



E 



/ n _i [u d (k+i) d \ £ fe _, ;/(fc) \ 

Y Y Y 1 ^{t)=y n 



Y J du d k d - 1 G{\\y n \\) + 0( 



n 



k=0 



u d n d G{a) + 0{n 2 



The calculation for the upper bound yields the same result. 

To finish the proof, we just need to bound the variance of these two 



random variables, using lemma 1372 



Var u n 



n-l M fe +!) 2 J 

Y Y Var (*i 

fc=0 j=\u> d k 2 \ 



■m(k),n\ n 



n-l [^(k+l?\ 

When d > 3, naming K d a suited constant depending only on the 
dimension d, we get: 

n— 1 



Var( Mn (^))~ < ^KJz^GiWvnW) 
Hence we get: 



k=0 



n d n d 



When d = 2, taking a suited constant C 2 gives us: 

n-l 

Var( Mn (^)) < ^«G(||^||)ln(n) 2 , 

k=0 

which yields: 



Va y[ u n (z)- \ < ln(n) 2 ^G(a) 



n 2 / n 2 



where C' 2 is a suited constant. 
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In both cases, the sum ^Var 



is finite. Since the same cal- 



culation can be applied to «+, we can use the same argument as in the 
proof of Theorem 12.11 to say that almost surely, as n tends to infinity, 



Mr. 



11" 



u d G(a). 



The proof of the second assumption is similar to that of the first one 
and relies on the following estimate: for d = 2, 

\ 



E 



Y Y Y 1 S i {t)=Vr. 

^k=0 j=lio 2 k 2 \ t=l 



J 



^4Hn 

fc=0 

An 2 In 



\Vn\ 



Id- 2 



\Vn\ 



k - 777(A)) 
+ 



+ 



+ 



\\Vn\\ n 



WnW n 



E 



And for d > 3, 

n-l L w<i ( fc+1 ) d J Zk-njik) \ 

Y Y Y 1 s i (t)=y n 



Y ^-^ d ~ x 



k=0 



d-2 



\Vn\ 



I d-2 



(k - r) I (k)) d - 



2n 2 
d-2 



?? 



\\Vr, 



d-2 



+ 

2 



n 



d-l 



ld-1 



O 



n 



d-l 



\Vn 



ld-1 



+ 



n 



The computations of the variances yield, for d = 2, 

n 



Var \ u n ( — J J < J2n 2 (Inn) 2 In ( n 
v Vn// WlVnW 



and for d > 3, 



Var ( u 



n 



< J d n 2 
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2-d 



where J2 and Jd are suitable constants depending only on d. It follows that 



Var 



Var 



M£)) 



n 2 In 



IWnll 



< 



n 2 In 



llVhll 



n 2 In 



llVnll 



< Id (\\Vr, 



\ n 



d-2 



These quantities have finite sums, which proves our result. 
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5 Back to a result of Lawler, Bramson and 
Griffeath 



The result in Theorem 12.11 can be used to derive a different proof of an 
existing time-scale result by Lawler, Bramson and Griffeath. 

In [2], the authors are interested in the time it takes to build a cluster 
of radius n, that is to say the total number of steps done by random walks 
during the construction of the cluster. They define A(t) as the internal 
DLA cluster on the time scale of individual random walks. Namely, if 

x{t) = max{/c such that <j\ + • • • + < t} 

then A(t) = A( X {t)). 

Theorem 5.1 (Lawler, Bramson, Griffeath, 1992) Let A(i) be the in- 
ternal DLA cluster on the time scale of individual random walks, then for 
all e > 0, almost surely, 

n(l+e)> 

where 

tn = n d+ * dUd 



d + 2 



While the original proof of this theorem studies the time spent by a 
given random walk inside the cluster, we have studied the time spent in 
a given point by all the random walks. This leads naturally to a different 
approach of the problem, in which we will sum the odometer function over 
all the points in the cluster. 

Let us call t' n the time taken to build the cluster A(n). Then 

t 'n= Yl Un ^ 

z:nzeZ d nA(n) 

Then for all n, 

U n(z) <f n < U ^ Z ) 

We can compute the following asymptotics, as n tends to infinity: 

o{n d+2 ) 



Y u n {z) = Y duJ dk d ' 1 u n i-J + 

eB„_„,^ k=i ^ ' 

k=i \ ^ 



z, nz eB n - Vl ( n ) 



d tu d n 2 > k f - + o(n a+ ), almost surely. 
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where / is defined by (12 .21) . Hence: 




u n (z) = d u d ri 



d+2 



I 



x d - 1 f(x)dx + o(n d+2 ) 



u n (z), so that we have, 
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